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WIMPs with electroweak scale masses (neutralinos, etc.) remain in kinetic equilibrium with other particle 
species until temperatures approximately in the range of 10 MeV to 1 GeV, leading to the formation of dark 
matter substructure with masses as small as 10~ 4 Mq to 10 -12 Mq. However, if dark matter consists of 
particles with MeV scale masses, as motivated by the observation of 51 1 keV emission from the Galactic Bulge, 
such particles are naturally expected to remain in kinetic equilibrium with the cosmic neutrino background until 
considerably later times. This would lead to a strong suppression of small scale structure with masses below 
about 10 7 Mq to 10 4 Mq . This cutoff scale has important implications for present and future searches for faint 
Local Group satellite galaxies and for the missing satellites problem. 

PACS numbers: 95.35.+d;95.30.Cq; FERMILAB-PUB-07-064-A; MADPH-07-1482 



Introduction. In the standard cosmology, featuring cold, 
collisionless dark matter, structures form hierarchically with 
the smallest mass objects forming first and progressively 
larger objects forming via subsequent mergers and accretion. 
This paradigm has been remarkably successful at describing 
the observed large scale structure. 

Within this paradigm, the mass of the smallest dark matter 
halos depends on the mass of the dark matter particles, and 
on the temperature at which they decouple kinetically from 
other particle species. For a typical weakly interacting mas- 
sive particle (WIMP) with an electroweak scale mass, kinetic 
decoupling from the standard model leptons occurs at a tem- 
perature in the range of roughly 10 MeV to 1 GeV, leading to 
the formation of structures with masses as small as 10" 4 M 
to 10~ 12 Mq If dark matter consists of particles which 

remain in kinetic equilibrium with neutrinos until later times 
(lower temperatures), the smallest dark matter halos will be 
considerably more massive than are predicted for WIMPs with 
electroweak scale masses J3]. 

Dark matter particles with MeV scale masses have been 
previously motivated by the observation of 511 keV emis- 
sion from the Galactic Bulge In particular, annihilat- 
ing MeV dark matter can inject the required rate of positrons 
into the Galactic Bulge, and also be produced in the early uni- 
verse with the measured dark matter abundance [5] . Dark mat- 
ter in the form of MeV mass scalars, <j>, annihilating through 
the exchange of a light gauge boson, U, can accommodate 
these requirements 01. Constraints on this scenario have 
been placed by colliders laM, neutrino experiments jlfJl . 
atomic phy sics experiments 111 ill , observations of supernova 
1987A 112], the 511 keV line width Q and Big Bang Nu- 
cleosynthesis (BBN) iQ. 

In this letter, we revisit the MeV dark matter scenario, and 
calculate the resulting matter power spectrum. We find that 
if the [/-boson's couplings to neutrinos is similar to its cou- 
plings to electrons, the matter power spectrum is suppressed 
on small scales, leading to an absence of dark matter halos 
with masses below about 10 7 Mm to 10 4 Mm. We show that 



this suppression scale is consistent with the region of param- 
eter space where the MeV dark matter has the correct relic 
abundance as determined by recent cosmological observations 

iB. 

The MeV Dark Matter Power Spectrum. To calculate the 
power spectrum for MeV dark matter, we start by determin- 
ing the temperature at which kinetic decoupling occurs. The 
squared amplitude for dark matter-neutrino elastic scattering 
is given by 
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where t = —2E^(1 — cos 9), m$ is the mass of the dark mat- 
ter particle, mu is the mass of the exchanged boson, and gu<t><t> 
and g\j vv are that boson's couplings to dark matter and neu- 
trinos, respectively. This leads to an elastic scattering cross 
section (in the mjj ^> E v limit) of 



0-4,1, 



„2 _2 p2 

yu<t><f>yuuv J2j L' 



(2) 



To determine the temperature of kinetic decoupling for the 
dark matter particle, we solve the Boltzmann equation, includ- 
ing the <j) — v collision term. The resulting equation lllrjll is 
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where F = 31n 3 gl <M> 9uvv T ?/(* 2m u m <t>) is the rate for 
the dark matter distribution function to relax to its equilib- 
rium value. An intuitive approximation to this relaxation rate 
is E k /E k = (An)- 1 JdE v dl{dri v /dE v )(dT^,/dl)SE k /E k 
where dn v jdE v is the differential number density of all neu- 
trinos, Ek = \p\ 2 /2m<f> is the kinetic energy of dark matter 
particles and SE k is the kinetic energy transferred per colli- 
sion. This approximation yields lb^g^^g^^T^ /(m^) 
for Ek = 3T„/2, which is a factor of about 1.5 smaller the 
exact result. 

Eq.[3]is solved by a Boltzmann distribution with a temper- 
ature for the dark matter particle that scales as T v during the 
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strongly coupled regime and as T^ after decoupling. We de- 
fine the kinetic decoupling temperature as T k d = T v , such 
that T{T v ) = H(T v )k 5.97 y/U^T^. This gives us 



T kd = 2.1 keV 
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Using the formalism of Ref. 101, we can use this result to cal- 
culate the power spectrum of MeV dark matter. 

The interactions and the subsequent decoupling of the dark 
matter particles leads to the damping of the matter power spec- 
trum. This results from three distinct processes. First, the 
coupling of the dark matter to other particle species intro- 
duces damped oscillatory features yfl. This is the dominant 
effect for the case of WIMPs with electroweak scale masses. 
Second, after decoupling, the free-streaming of the dark mat- 
ter particles further suppresses the power spectrum. For MeV 
dark matter, this effect dominates for the viable region of pa- 
rameter space where T^d > keV. Third, as neutrinos kineti- 
cally decouple from the dark matter they begin to free-stream 
and damp the power spectrum further. This effect, however, is 
subdominant. 

In Fig. Q] we show the effect on the matter power spec- 
trum of MeV dark matter as compared to that for the stan- 
dard cold dark matter case. Large wavenumbers are strongly 
suppressed, resulting in reduced number of small dark mat- 
ter halos. Also shown in the figure as a dotted curve is the 
(strictest) limit found for the case of warm dark matter from 
observations of the lyman-alpha forest II 1711 . 

For Tkd <L keV, the scale at which the power spectrum is 
truncated is closely related to the free-streaming scale, 



,-i oci AeV\ 1/2 /McVV , 
k f = 2.5 kpc ( jtJ \~^) ln ( 4a E Q /akd) ; 
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where a^d and oeq are the scale factors at decoupling and 
matter-radiation equality, respectively. The suppression of 
the dark matter power spectrum on scales smaller than kj 1 , 
in turn, leads to a cutoff in the mass function of dark 
matter halos. Compared to the case with no cutoff, one 
would find a paucity of halos with masses less than roughly 
4ir(ir/kf) 3 pu /3, where pm is the present cosmological mat- 
ter density. To obtain a more accurate estimate, we find the 
mass at which the expected number of dark matter halos falls 
by a factor of e compared to the prediction for dark matter par- 
ticles with electroweak scale masses. We calculate the mass 
function of dark matter halos using the Press-Schechter pre- 
scription. We note that the validity of this prescription for 
power spectra with sharply truncated power (as found in our 
scenario) has not been conclusively demonstrated. Neverthe- 
less, the cutoff mass derived here is useful in the sense that 
it highlights the mass scale below which we expect deviations 
from the predictions of standard cold dark matter. We find this 
cutoff mass to be 
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FIG. 1: The effect of keV scale kinetic decoupling on the matter 
power spectrum, as predicted in MeV dark matter. Shown are results 
for a 1 MeV dark matter particle with a 10 keV (solid) and 1.0 keV 
(dashed) kinetic decoupling temperature. The dotted line denotes the 
limit relevant for warm dark matter, as inferred from observations of 
the lyman-alpha forest fTvh . 



Combining this expression with our specific particle physics 
scenario, we arrive at the estimate: 



M c ~ 10 7 M P 
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We note that for Tkd ~ keV and ~ MeV, the small- 
est halos that form (those with mass ~M C ) are the ones that 
host the smallest of the dwarf galaxies seen in the Milky Way 
1 18D . Therefore, the predictions for the number of Milky Way 
satellites will be different in this scenario compared to that 
for dark matter with electroweak scale masses. Numerical 
simulations with truncated power spectra that are able to re- 
solve halos with masses below M c and a detailed treatment 
of galaxy formation on small scales will be required to make 
robust predictions for the satellite (dwarf) galaxy population 
in galaxies like the Milky Way and Andromeda. 

Relic Abundance and Other Constraints. There are a 
number of constraints on the various couplings and masses in 
MeV dark matter. First, we require that dark matter is ther- 
mally produced in the early universe with the observed abun- 
dance [ 15]. The annihilation cross section for scalar dark mat- 
ter particles through the s-channel exchange of a [/-boson is 
Si: 
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FIG. 2: The thermal relic abundance of dark matter as a function of 
its mass for ray = 1 MeV (solid), 3 MeV (dashed) and 10 MeV 
(dotted). In each case, the product gu<p<i>guff was set to 10~ 6 for 
each of / = e, v e , z/ M and v T . The dashed horizontal lines denote 
the measured density of dark matter 1 15]. For illustration, we have 
shown the complete range for the dark matter mass, although only 
iricf, < rau is viable 1 1911. 



where we have denoted the U couplings to left and right- 
handed fermions by gf L and gf R , respectively. The sum is 
over e + e~ and the three species of neutrinos. Notice that at 
low velocities (s » 4mJ) the cross section approaches zero, 
being entirely the result of a P-wave amplitude. 

In order for dark matter annihilations to generate the ob- 
served 511 ke V photons from the Galactic Bulge, positrons 
must be injected with energies no greater than ~ 3 MeV (more 
energetic positrons would unacceptably broaden the 511 keV 
line width). This leads to the constraint, 0.511 MeV < < 
3 MeV [13]. We also require that < mu in order to avoid 
dark matter annihilating largely to UU, which is not s-wave 
suppressed 111911 . For a 0.511-3 MeV dark matter particle to 
be generated in a quantity consistent with the observed dark 
matter abundance, an annihilation cross section on the order 
of a picobarn is required during the freeze-out epoch. 

We show in Fig. [2] the abundance of dark matter in this 
model as a function of its mass, for three values of the gauge 
boson mass, and for couplings of gu^4> x 9uff = 10~ 6 . 
Throughout, we adopt a common [/-fermion-fermion cou- 
pling for electrons and neutrinos. 

The [/-boson's couplings to fermions are constrained by ve 
scattering experiments such that g Uvv \J g UeLEL + 9ue R e R ^ 
mfjGF 0,0] • For the case of a common [/-fermion-fermion 
coupling, this reduces to guff Si 2.9 x 10~ 6 x (my /MeV) 2 . 
A somewhat weaker constraint can be found from measure- 
ments of the electron's magnetic moment 0. 

In Fig. |3] we show the range of rau an d the product gu<t>4> x 
guff for which the measured dark matter density can be made 



to match the thermal relic abundance in this model (for some 
value of rricf) in the range of m e to 3 MeV). We also show the 
constraints from ve scattering experiments (for the optimal 
case of gucjicji ~ 1) 00] and from the measurement of the 
electron's magnetic moment 0. 

As light blue lines, we have plotted contours of constant M c 
from 10 7 to 10 4 solar masses. Here, we have used — 1 
MeV. For other masses, the results vary as M c cx m^ 15 ^ 8 . 
From this figure, we see that once all of the constraints are 
considered, M c is generally expected to fall in the range of 10 4 
to 10 7 Mq. It should be noted that the region where mu ~ 
1 — 6 MeV and gu<j><t>9uff Si 10~ 7 is highly fine tuned and 
relies on being very close to the resonance at 2m,f > ~ mu to 
avoid the overproduction of dark matter. 

The coupling of MeV dark matter to the neutrinos and elec- 
trons could change BBN predictions by causing the neutrino 
and photon temperatures to be the same to lower redshifts than 
is standard. Effects of this nature have been considered pre- 
viously 11411 . though it is unclear how these constraints apply 
to this model on account of the presence of extra thermalized 
scalars and new neutrino interactions; a detailed analysis of 
the constraints from BBN is beyond the scope of this letter, 
but we have checked that there are viable regions of param- 
eter space where the expected deviations from standard BBN 
predictions are within observational bounds J20ll . 

MeV dark matter (and associated [/-boson) with couplings 
to neutrinos would have other observable consequences. The 
existence of such a [/-boson would lead to TeV scale ab- 
sorption features in the high-energy cosmic neutrino spec- 
trum 112 If . The spectrum of neutrinos produced in core- 
collapse supernovae could also be modified due to their in- 
teractions with dark matter particles produced during the col- 
lapse Il2ll . 

Conclusions. In summary, we have calculated the small 
scale power spectrum of MeV dark matter. This scenario is 
motivated by the observation of 5 1 1 keV emission from the 
Galactic Bulge. Assuming that the relevant couplings to neu- 
trinos are similar to those to electrons, we find that MeV dark 
matter particles remain in kinetic equilibrium with the cosmic 
neutrino background up to temperatures of ~1-10 keV. This 
late kinetic decoupling leads to larger free-streaming lengths 
for MeV dark matter as compared to WIMPs with electroweak 
scale masses. This highly suppresses the formation of small 
scale structure. Depending on the parameters considered, the 
matter power spectrum is expected to be cutoff below 10 7 Mq 
to 10 4 Mq in this scenario. 

This result has a number of particularly interesting astro- 
physical implications. First, it predicts a cutoff in the mass 
function of dwarf galaxies at a mass scale much larger than 
that for WIMPs with electroweak scale masses. It was previ- 
ously shown that the number of dwarf galaxy-sized dark mat- 
ter halos in numerical simulations of cold dark matter is con- 
siderably larger than the observed populations in the Milky 
Way and Andromeda galaxies (ie. the "missing satellites prob- 
lem") B221 12311 . This issue may be resolved by astrophysical 
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FIG. 3: Regions in the mu versus gu<f»j>9uff plane in which the 
measured dark matter density matches the thermal relic abundance 
for some value of in the range of m e to 3 MeV. We have adopted 
a common U -fermion-fermion coupling for electrons and neutrinos. 
The dashed line denotes the constraint from ve scattering experi- 
ments (for the optimal case of gu<t>4> ~ 1) la |2D- The dotted line 
denotes the (weaker) constraint from measurements of the electron's 
magnetic moment [6]. The light blue lines are contours of constant 
M c from 10 4 to 10 7 solar masses. Here, we have used = 1 MeV. 
From this figure, we see that once all of the constraints are consid- 
ered, M c in the range of 10 4 Mq to 10 7 Mq are generally expected. 

means [24] or by altering the nature of the dark matter's inter- 
actions 125| or mechanism of production [26]. The model we 
consider falls into this latter category, though detailed numer- 
ical simulations will be required to make precise predictions 
of the satellite population in MeV dark matter. 

Tests of dark matter models with observations of small 
scale structure, as we have discussed in this paper, are becom- 
ing a more realistic possibility given the recent discoveries of 
faint satellite companions to the Milky Way and Andromeda 
galaxies 02711 . Present estimates of the masses of these new 
satellites fall in the range 10 5 M Q to 10 7 M , which is near 
the cutoff mass scale in MeV dark matter. Present and fu- 
ture searches for faint satellites, and the characterization of 
the mass function at these scales, will thus provide important 
constraints on MeV dark matter. 
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